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Abstract 


In  order  to  model  high-strain-rate  deformation  and  texture  evolution 
in  commercially  pure  tantalum  (Ta),  a  description  for  the  thermal- 
elastic-viscoplastic  behavior  of  Ta  single  crystals  is  considered  along 
with  an  associated  polycrystal  averaging  scheme.  The  description  incor¬ 
porates  a  temperature-dependent  model  for  pencil  glide  on  the  planes  of 
maximum-resolved  shear  stress.  Calculated  stress-strain  data  and  tex¬ 
ture  evolution  for  this  model  are  compared  to  those  of  a  restricted-glide 
model  and  to  experimental  data. 
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1.  Introduction 


Tantalum  (Ta)  is  a  high-density  body-centered  cubic  (bcc)  metal  that  is  highly  ductile 
over  a  wide  range  of  temperatures  and  strain  rates.  The  high  density  of  Ta  makes  it  an 
attractive  candidate  for  penetration  applications  while  its  high  ductility  translates  to  high 
formability  under  the  broad  range  of  processing  environments  that  must  be  tolerated  by  an 
explosively  deformed  liner.  Wrought  Ta  is  most  commonly  available  in  the  form  of  rolled 
plate,  which  must  undergo  deep  drawing  to  form  the  initial  liner.  Prior  to  the  penetration 
event,  the  initial  liner  is  either  deep  drawn  or  fully  collapsed  at  high  strain  rate  to  form  an 
elongated  cup  or  rod  in  the  reverse  direction  of  initial  drawing. 

The  drawability  of  sheet  material  has  long  been  recognized  to  depend  on  material  flow 
coming  from  the  width  and  not  the  thickness  of  the  rolled  sheet  (Backofen  1972).  For  bcc 
steels,  this  flow  property  has  been  shown  to  depend  on  the  alignment  of  {111}  planes  with 
the  rolling  plane  of  the  sheet  (Lankford,  Snyder,  and  Bauscher  1950).  Such  crystallographic 
texture  considerations  are  also  expected  to  influence  the  deep-drawing  rolled  Ta  sheet.  With 
this  in  mind,  Clark  et  al.  (1991)  showed  the  effects  of  forging  and  rolling  conditions  on 
through-thickness  texture  gradients  in  Ta  liner  material.  Their  goal  was  to  produce  Ta  sheet 
with  a  homogeneous  through-thickness  microstructure. 

In  order  to  develop  predictive  capabilities  for  the  thermomechanical  processing  of  Ta,  a 
temperature-  and  rate-dependent  model  is  presented  for  flnite-strain  deformation  by  <  111  > 
crystallographic  slip  that  is  not  restricted  to  particular  crystallographic  planes  (pencil  glide). 
The  model  is  based  on  the  Asaro  and  Needleman  (1985)  elastic-viscoplastic  formulation  for 
rate-dependent  crystallographic  slip  in  polycrystalline  aggregates.  In  this  model,  the  stress 
state  in  each  particular  grain  is  uniquely  determined  from  the  elastically  deforming  crystal 
lattice.  This  stress  state  is  then  used  to  determine  crystallographic  planes  via  a  maximum- 
resolved  shear-stress  criteria  originally  outlined  by  Becker  (1995).  The  determination  of 
the  slip  planes  is  conducted  throughout  the  entire  deformation  process  (i.e.,  evaluated  at 
each  time  step  during  the  calculation).  This  model  allows  for  the  wavy  glide  nature  of 
<  111  >  {110}  and  <  111  >  {121}  crystallographic  slip  in  bcc  metals.  It  should  also  be 
noted  that  although  only  homogeneous  deformations  axe  considered  at  present,  this  model 
may  be  particularly  well  suited  for  implementation  within  a  flnite-element  framework  to 
simulate  the  entire  thermomechanical  process.  Since  the  maximum-resolved  shear-stress 
criteria  determines  the  proper  slip  planes  as  the  computation  proceeds,  there  is  only  a  single¬ 
plane  normal  calculated  per  slip  direction,  and  the  orientation  of  this  normal  is  allowed  to 
evolve  during  the  deformation.  This  approach  reduces  the  number  of  slip  systems  typically 
required  to  model  bcc  restricted  glide  to  only  the  four  pencil-glide  directions  and  their 
associated  slip-plane  normals.  The  temperature  dependence  of  the  slip  rate  is  incorporated 
by  the  exponential  temperature  dependence  recently  proposed  by  Lee  et  al.  (1997).  As  such, 
the  model  can  account  for  thermal  softening  due  to  adiabatic  heating  that  accompanies 
high-rate  deformation  during  explosive  forming.  The  stress-strain  behavior  and  associated 
texture  evolution  predicted  by  the  model  are  compared  with  experimental  results  obtained 
at  high  strain  rates. 

Notation  used  in  this  paper  is  based  on  the  following  conventions.  Scalars  are  denoted  by 
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mathematical  italics  or  calligraphics  {A,  a,  a,  X),  vectors  by  lowercase  boldface  (a),  second- 
order  tensors  by  uppercase  boldface  (A),  and  fourth-order  tensors  by  calligraphic  uppercase 
boldface  (K).  The  boldface  Greek  letters  designate  second-order  tensors.  The  summation 
convention  is  used  for  Latin  indices.  Superscripts  pertaining  to  slip  systems  are  written  in 
Greek  letters;  summation  over  these  are  indicated  explicitly.  Single  dots  and  double  dots  are 
used  to  indicate  the  following  products,  respectively: 

F  •  m  ^  Fijrrij; 

/C  :  L  ^  JCijkiLik  &nd  t  •  D  ^  Tik^kj  • 

The  dyadic  product  of  two  side-by-side  tensors  yields  a  tensor  with  a  rank  equal  to  the  sum 
of  the  ranks;  an  example  is  sm  ^  Si-rrij.  Finally,  time  derivatives  are  denoted  by  superposed 

dots,  e.g., 

F  =  dF/dt . 

2.  Constitutive  Relations  for  Single  Crystals 


The  description  of  the  rate-dependent  finite-strain  behavior  of  single  crystals  deforming 
by  crystallographic  slip  accompanied  by  elastic-lattice  distortion  presented  here  follows  from 
the  work  of  Asaro  and  Needleman  (1985). 


2.1  Single- Crystal  Kinematics 

Plastic  deformation  is  assumed  to  occur  by  simple  shears  along  the  various  slip  systems  of 
the  crystal.  This  shear  is  described  in  a  bulk  continuum  framework  and  is  due  to  dislocation 
movement  along  the  slip  systems  or  other  mechanisms  that  produce  simple  shear.  In  addition 
to  the  plastic  shearing,  the  lattice  undergoes  elastic  distortion  and  rigid-body  rotation.  The 
motion  of  the  single  crystal  can  then  be  decomposed  into  a  plastic  part,  which  de^ribes 
crystallographic  glide,  and  an  elastic  part,  which  contains  lattice  stretching  and  rotation.  If 
F  designates  the  deformation  gradient  of  the  single  crystal,  a  multiplicative  decomposition 
of  the  type  due  to  Lee  (1969)  is  written  as 


F  =  F*FP, 


(1) 


where  F^  is  due  to  the  net  motion  of  dislocations,  and  F*  contains  the  elastic  distortion  and 
rigid-body  rotation  of  the  lattice  (Rice  1971;  Hill  and  Havner  1982;  Asaro  1983).  The  slip 
direction  of  a  particular  system,  a,  is  denoted  by  and  must  lie  in  the  slip  plane  with 
normal  The  orthonormal  pair  and  convect  with  the  lattice  to  become 


g*(Q:)  _  F*.s(“) 


(2) 


Note  that  the  total  velocity  gradient  L,  can  be  decomposed  into  its  elastic  and  plastic  parts, 
L*  and  L^,  as 

L  =  F  •  F"^  =  L*  -h  LP  .  (3) 
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The  plastic  part  of  the  velocity  gradient  can  be  decomposed  into  its  symmetric  and  skew 
symmetric  parts  and  respectively,  as 

LP  =  DP  +  .  (4) 


We  incorporate  equation  (2)  into  the  single-crystal  kinematics  by  writing  Lp  in  terms  of  the 
shear  rates  on  each  system.  Thus, 


a 

(5) 

where  a  is  summed  over  all  active  slip  systems,  and  7^“^  is  the  rate  of  shearing  on  slip  system 
a.  For  each  slip  system,  the  symmetric  and  antisymmetric  parts  of  the  Schmid  tensor  are 
defined  as 

p(“)  =  (6) 

2 

and 

2 

(7) 

so  that 

DP  = 

a 

(8) 

and 

(9) 

a 


The  total  stretch  tensor  D  and  the  total  spin  are  decomposed  as  D  =  Dp  -f  D*  and  fZ  = 
fiP  4-  D*,  where  D*  is  the  elastic  stretch  tensor,  and  D*  is  the  lattice  spin. 

2.2  Crystal  Elasticity 

V 

If  we  denote  the  Kirchhoff  stress  as  r,  the  Jaumann  rate  of  this  stress,  t,  on  axes  that 
spin  with  the  material  is  ^ 

T  =  T—  rt’T+T’Ct.  (10) 

V 

We  can  also  write  the  Jaumann  rate  of  Kirchoff  stress  ,  t*,  formed  on  axes  that  spin  with 
the  lattice  as 

T*  =  -i— n*-T-f-T-D*.  (11) 

As  shown  by  Asaro  and  Needleman  (1985),  for  a  material  characterized  by  an  elastic-strain 
energy  function  (elastic  moduli,  £),  r*  is  given  by 

t*  =  C:D*  +  D*t+tD*^C:B*.  (12) 

This  approximation  is  due  to  the  assumption  that  stresses  are  small  compared  to  elastic 

moduli.  We  incorporate  equations  (9)  and  (10)  into  equation  (11)  to  get 

r*=T  +  •  T - T  •  .  (13) 
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Equation  (13)  is  now  substituted  into  equation  (12).  The  resulting  constitutive  law  is  written 
as 

T  =C  :D  -  + 

a 

W(°)  .  (14) 

This  expression  is  often  written  in  more  concise  form  as 

T  =  r:D-i:R(“>7^“^  (15) 


2.3  Slip  Constitutive  Model 

To  capture  the  dynamic  behavior  of  Ta  at  high  strain  rates,  the  slip  constitutive  law 
at  time  t  will  depend  on  the  current  stress  state  and  system  resistance  to  glide,  t  and 
respectively,  as  well  as  the  current  temperature,  T,  in  the  deforming  crystal.  Generally,  the 
shear  rate  7^“^  of  the  system  a  can  be  specified  by  a  function  with  the  form 

^(“)(t)  =  f“^(^,5(“),T).  (16) 

The  specific  form  that  we  will  examine  is  a  modification  of  a  widely  used  power  law  due  to 
Lee  et  al.  (1997)  that  relates  these  shear  rates  to  their  conjugate  measure  of  stresses  (i.e., 
resolved  shear  stresses)  viz., 

v(o)  \ 

—  I  5^71(7^“^)  exp(—AAr),  (17) 

To  / 

where  70  is  a  reference  shear  rate  assumed  the  same  for  all  slip  systems,  m  is  the  strain- 
rate  sensitivity  coefficient,  A  is  the  thermal  softening  coefficient,  and  AT  is  the  temperature 
increase  from  a  reference  temperature  (T  -  T^e/).  The  reference  temperature  for  this  work 
is  taken  as  ambient  room  temperature.  Initially,  each  g^°‘^  is  taken  to  be  a  constant  that  can 
be  different  from  one  slip  system  to  another,  and  the  rate  of  increase  of  these  resistances  is 
specified  by  the  following  hardening  law: 

^(“)  =  ^h“^|7(^)|,  (18) 

0 

and 

h^0  =  (no  sum  on  j3),  (19) 

where  is  a  matrix  fully  populated  by  ones.  The  single-slip  hardening  rate,  h^,  will  be 
given  the  saturation  form  proposed  by  Kalidindi,  Bonchorst,  and  Anand  (1992); 

=  (20) 
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where  g'^  is  the  saturation  resistance  to  slip.  The  expression  for  the  dissipation  of  mechanical 
work  into  heat  during  high-strain-rate  deformation  is  given  as 


AT= 

•/o  pO.  V 


(21) 


where  0  <  x  ^  1  represents  the  fraction  of  plastic  work  that  is  converted  to  heat,  p  is  the 
density,  and  Cp  is  the  heat  capacity. 


2.4  Determination  of  the  Slip  Planes 

In  this  work,  high  rate  deformation  in  Ta  is  assumed  to  occur  by  pencil  glide  along  the 
<  111  >  crystallographic  directions.  The  slip  planes  will  not  be  restricted  to  any  specific 
set  of  crystallographic  planes,  but  will  be  allowed  to  occur  on  any  plane  that  maximizes  the 
resolved  shear  stress  for  a  particular  slip  direction.  The  following  method  has  been  proposed 
by  Becker  (1995)  for  the  determination  of  the  plane  that  maximizes  the  resolved  shear  stress. 
For  a  particular  slip  system,  a,  the  resolved  shear  stress  t°  is  given  as  a  function  of  the  slip 
geometry  and  the  applied  stress 

^.(q)  _  s*(a)  _  .  ^*(q)  ^22) 

The  vector  represents  the  traction  normal  to  the  slip  direction  in  the  updated  configu¬ 
ration. 

Note  that  r®  is  maximum  when  lies  in  the  direction  of  The  slip-plane 

normal,  therefore,  is  calculated  by  removing  the  component  of  in  the  direction  of  s*^®^ 
(to  assure  orthogonality  between  s*^®)  and  m*^®^).  Normalization  of  the  the  slip-plane  normal 
then  yields 

The  plane  normal  is  computed  for  each  of  the  four  slip  directions.  Thus,  only  four  slip 
systems  are  selected  instead  of  12,  24,  or  48,  as  in  the  case  of  restricted  glide;  this  results  in 
an  important  computational  time  reduction. 


2.5  Numerical  Implementation. 

During  the  computations,  care  must  be  taken  in  order  to  estimate  the  actual  slip  rates 
over  a  given  time  increment.  In  order  to  obtain  a  conditionally  stable,  explicit  integration 
of  the  stress  rate  across  a  time  step,  the  tangent  modulus  method  of  Peirce,  Asaro,  and 
Needleman  (1983)  is  employed.  The  following  is  a  modification  of  their  method  in  order 
to  accommodate  the  temperature-dependent  slip-rate  constitutive  model.  By  expanding 
equation  (16)  in  a  Taylor  series  about  the  current  time  t,  we  can  estimate  the  increment  in 
slip  across  the  current  time  step. 

7®(i  -I-  At)  =  7®(t)  4-  -^"Ag®  + 


5 


(24) 


We  insert  equation  (17)  into  equation  (24)  to  yield 

'y^{t  +  At)  = 

[l  +  ^(^  +  1f +-'^T)]7«(i).  (25) 

The  increment  of  slip  across  the  time  step,  defined  by 

A7“  =  7“(t  +  At)-7“(<),  (26) 

is  approximated  by  Pierce,  Asaro,  and  Needleman  (1983)  as  a  linear  interpolation  across  the 
time  step 

A7“  =  [(1  -  0)r{t)  +  Or{t  +  At)]  At,  (27) 

where  0  is  an  interpolation  parameter  taken  here  as  0.5.  Substitution  of  equation  (25)  into 
the  previous  equation,  equation  (27),  yields 


A7“  = 


7“(i)  + 


,  Ag‘ 


m 


4-  +  AAT 


At . 


(28) 


Using  At  =  fAt,  Ag  =  gAt  and  equation  (21)  for  AT  in  equations  (28)  and  (27)  we  get 


A7“  =  7“(t)At  + 


eAtY^it) 


mr^ 


R“  :  D  -  S  (r-  :  P'  +  ^(t) 


(29) 


This  estimate  for  the  increment  of  slip  across  the  time  step  is  now  simplified  and  made 
consistent  with  the  original  Pierce,  Asaxo,  and  Needleman  (1983)  notation.  By  expressing 


/eAtrjt) 

\  rriT^ 


Q"  = 


R“ 


(30) 


and 


Nafi  =  5q/3  + 


y  mT° 


(O' 


R“  :  Axr^ 


9^  pCv  J 


(31) 


We  can  now  substitute  equation  (29)  into  equation  (15)  for  a  stable  estimate  of  the  stress 


rates  across  the  time  step; 

with 


r  =  C:D-5;R“r, 

a 

c  =  c-Y, 


(32) 


=  Y 


r  =  Y  . 


(33) 
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The  matrix  M  is  the  inverse  of  N  (guaranteed  invertible  for  sufficiently  small  time  steps) 
(Pierce,  Asaro,  and  Needleman  1983). 

2.6  Polycrystal  Averaging 

The  modified  Taylor  polycrystal  averaging  scheme  (e.g.,  Asaro  and  Needleman  1985) 
is  incorporated  in  order  to  (1)  localize  the  macroscopic  deformation  to  the  local  grains 
and  (2)  average  the  resulting  stresses  from  the  local  grains  to  the  associated  macroscopic 
aggregate  stresses.  The  model  assumes  that  deformation  within  a  grain  is  uniform  and 
that  each  grain  undergoes  the  same  deformation  as  the  entire  aggregate.  Therefore,  the 
local  velocity  gradient,  L,  is  identical  to  the  macroscopic  velocity  gradient,  L.  In  this  way, 
deformation  compatibility  between  the  grains  and  macroscopic  deformation  is  maintained. 
After  localization  of  the  applied  L  to  L,  the  local  Cauchy  stresses,  S,  are  then  calculated 
and  volume-averaged  to  achieve  the  macroscopic  stress,  S.  Since  the  calculations  assume 
that  all  the  grains  occupy  an  equal  volume  within  the  polycrystal, 

s  =  if;s'»>,  (34) 

9=1 

where  K  is  the  number  of  grains  in  the  aggregate. 


3.  Results  and  Discussion 


Lee  et  al.  (1997)  compressed  samples  of  pure  Ta  at  strain  rates  ranging  from  0.001/s  to 
30, 000/s  and  temperatures  ranging  from  25°C  to  525"C.  The  previously  described  model  for 
pencil- glide  has  been  calibrated  to  their  room- temperature  data  at  a  strain  rate  of  3, 000/s  via 
the  folloAving  procedure.  The  elastic  constants  were  taken  from  Gray  and  Vecchio  (1995)  as 
Cn  =  260.9  GPa,  C12  =  15.74  GPa,  and  C44  =  81.8  Gpa.  The  thermal-softening  parameter, 
A,  was  set  to  0.00153A~^  to  reflect  the  thermal  softening  observed  from  Lee  et  al.  (1997) 
between  their  quasi-static  room-temperature  data  and  data  taken  at  325"C.  The  portion 
of  mechanical  work  converted  to  heat,  X)  was  assumed  to  be  0.9.  The  initial  resistance  to 
slip,  5o  was  estimated  at  154  MPa  and  taken  to  be  equal  on  all  slip  systems,  while  the 
reference  shear  rate  was  taken  to  be  70  =  1000/s.  The  strain-rate  hardening  exponent  was 
taken  as  m  =  0.084  accounting  for  the  strain-rate  sensitivity  between  their  quasi-static  and 
3,000/s  data  at  an  offset  strain  of  0.05.  Since  the  strain-hardening  behavior  at  large  strains 
(>  0.2)  is  quite  different  between  these  two  sets  of  data,  it  is  not  expected  that  a  single  set 
of  parameters  will  model  this  entire  data  set  for  the  current  model.  The  remaining  set  of 
material  parameters  were  fit  using  the  minimization  routine  “e04upf  ’  in  the  NAG  library  in 
order  to  minimize 


where 


$(x)  =  _f(x) .  f(x), 

1  ^  {^22)  _ 

f(x)  ^  i^experiment  ^ 


7 


over  n  data  points. 

The  components  of  x  are  taken  as  the  strain-hardening  constants  for  the  slip  constitutive 
law;  ho,  gf,  and  a.  ^22 (£22)  is  the  component  of  stress  along  the  compression  direction  for 
the  sample.  Figure  1  shows  the  resulting  fit  for  uniaxial  compression  at  3,000/s.  The 
strain-hardening  constants  are  given  ss  ho  =  2498.6  MPa,  gf  =  231.6  MPa,  and  a  ==  1.08. 

The  deformation  texture  predicted  from  this  calibrated  model  is  compared  to  equivalent 
textures  resulting  from  restricted-glide  on  the  <  111  >  {110}  and  <  111  >  {121}  systems 
in  Figure  2.  After  a  true  strain  of  0.70,  both  models  compare  well  to  the  experimentally 
measured  textures  found  in  Lee  et  al.  (1997)  and  are  nearly  indistinguishable  from  one 
another.  It  is  important  to  note  here  that  the  restricted  glide  calculation  over  24  shp  systems 
required  almost  20  times  more  computational  effort  than  the  pencil-glide  calculation  over 
four  slip  systems;  yet,  both  produce  macroscopically  equivalent  grain  distributions. 


Figure  1.  Pencil-glide  model  (solid  line)  calibrated  to  the  stress-strain  data  of 
Lee  et  al.  (1997)  (dashed  line). 

In  order  to  understand  how  these  equivalent  textures  evolve,  a  single  grain  was  simulated 
to  imdergo  uniaxial  compression  along  the  crystallographic  <  100  >,  <  110  >  and  <  111  > 
directions.  Table  1  summarizes  the  calculated  slip-plane  normals  at  the  onset  of  plastic 
straining  (a  Mises  effective  stain  of  «  0.005).  It  is  interesting  to  note  that  the  slip-plane 
normals  for  compression  along  the  <  100  >  and  <  110  >  directions  are  {110}  and  {121}, 
consistent  with  those  observed  in  bcc  materials.  During  compression  along  the  <  111  >,  the 
calculated  normals  are  very  close  to  the  {123}  planes,  which  are  also  observed  sfip  planes 
in  bcc  materials.  Since  the  restricted-glide  model  incorporated  the  <  111  >  {110}  and 
<  111  >  {121}  systems  explicitly,  calculated  grain  rotations  were  almost  identical  for  the 
two  models.  Perhaps  more  importantly,  however,  is  the  fact  that  the  pencil  glide  model  seeks 
to  maximize  the  resolved  shear  stress  along  the  slip  direction  based  on  the  overall  stress  state 
in  the  grain.  The  stress  response  in  the  grain  is  due  to  the  incorporation  of  the  single-crystal 
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Pencil 

glide 


Restricted 

glide 


(100)  Pole  figures  (110)  Pole  figures 

Figure  2.  Calculated  pole  figures  for  uniaxial  compression  of  300  initially 

random-oriented  grains  using  (a)  the  pencil-glide  model  with  four 

<  111  >  directions  and  (b)  a  restricted-glide  model  with  24 

<  111  >  {110}  and  <  111  >  {121}  systems. 

elasticity  and  is  therefore  perhaps  more  physically  meaningful  than  restricting  slip  to  only 
the  “observed”  slip  systems. 

In  closing,  it  should  be  noted  that  a  temperature-dependent  elastic- viscoplastic  model  for 
pencil  glide  has  been  introduced  and  shown  capable  of  representing  the  mechanical  behavior 
and  texture  evolution  in  pure  polycrystalline  Ta.  This  new  model  calculates  slip  planes 
so  as  to  maximize  the  resolved  shear  stress  in  the  slip  direction  via  the  stresses  due  to 
elastic-lattice  distortion  and  therefore  provides  a  physically  based  description  of  the  active 
deformation  mechanisms. 
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Compression 

Axis 

Slip 

Direction 

Plane 

Normal 

1  1  1 

1.00-2.00  1.00 

10  0 

-111 

-1.00-2.00  1.00 

-1-1  1 

-1.00  2.00  1.00 

1-11 

1.00  2.00  1.00 

1  1  1 

1.00-1.00  2.00 

1  1  0 

-111 

-1.00-1.00  0.00 

-1-1  1 

1.00  1.00  2.00 

1-1  1 

1.00  1.00  0.00 

1  1  1 

-2.00  1.00  1.00 

1  1  1 

-1  1  1 

-2.00  -1.00  -1.00 

-1-1  1 

1.78  1.00  2.78 

1-1  1 

-1.78  -2.78  -1.00 

Table  1.  Calculated  slip-plane  normals  as  a  function  of  compression  direction 
for  a  single  crystal  under  uniaxial  compression. 
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